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The retarded electron propagator Sn(po, p) at high momentum p> T was shown by Blaizot and 
Iancu to be an entire function of complex po- In this paper a specific form for Sr(po, p) is obtained 
and checked by showing that its temporal Fourier transform Sii(t,p) has the correct behavior at 
large t. Potential infrared and collinear divergences from the emission of soft photons do not occur. 
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I. INTRODUCTION 

High temperature QCD is, in one respect, simpler 
than high temperature QED. In both theories the damp- 
ing rates for charged particles at very high momentum, 
p 3> T, have infrared singularities at the mass shell 
when computed in perturbation theory using Braaten- 
Pisarski resummation |l| . The divergence comes from the 
emission and absorption of quasistatic, transverse gauge 
bosons. In QCD this apparent infrared divergence in the 
quark and gluon damping rates is cut off by the nonper- 
turbative magnetic screening mass for transverse gluons 
and this gives finite damping rates 0- For quarks the 
damping rate 7 determines the location of the pole in 
the retarded propagator sa(po,p): 



quark 



(po,p) 



p a - p + i-f/2 



The Fourier transform gives the time-dependent retarded 
propagator 



quark 
9 R 



(t,p) = -i6(t)e~ zpt exp[-7*/2] 



In QED there is no magnetic screening mass and con- 
sequently the lowest order electron damping rate is in- 
frared divergent 0, HH 0> IE E3 • A direct analysis of the 
time dependent retarded propagator sii(t,p) produced 
surprising results: the time dependent propagator is not 
infrared divergent and at large times, uj p t 3> 1, it falls 
more rapidly than an exponential 



sn(t,p) 



-ie lp *exp 



{-(m[ln(u p t) + C]\. (1.1) 



Here a = g 2 /Air is the fine structure constant and uj p — 
gT/3 is the plasma frequency. (It is customary to use g 
for the electric charge so as not to confuse it with the base 
of the natural logarithm e.) Blaizot and Iancu @)H,@ ob- 
tained this result using a functional integral version of the 
Bloch-Nordsieck approximation. This resums an infinite 
class of Feynman graphs with no fermion loops, except for 
the Braaten-Pisarksi hard-thermal-loop correction to the 
quasistatic transverse photon propagator. Wang, Boy- 
anovsky, de Vega, and Lee obtained the same result by 
using the dynamical renormalization group to improve 
the lowest order result [10( . Boyanovsky et at |llj found 
the same behavior for the retarded propagator of the 
scalar particle in scalar QED. Because su(t,p) falls faster 



than an exponential, its Fourier transform sr{j>q,p) has 
some very unusual properties. Blaizot and Iancu proved 
the following about sr(pq,p): (1) It is an entire func- 
tion of pa; (2) It vanishes in the limit Impo ~~ * +00; 
(3) It diverges in the limit Impo ~~ > —00. Specifically, at 
Po = p — i£ it is pure imaginary and bounded from below: 



1 



aT 



(—>oo: \s R (p-i(,,p)\ > — explaexp 

(1.2) 

This asymptotic behavior was deduced in |8| using C = 
and the constant a was shown to have the value a = 
3<?exp(— 1)/Att. However, revising the argument of @ 
for a non-zero value of C gives 
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a = — exp[— C — 1] 
4tv 



(1.3) 



A propagator that is an entire function of po has no poles, 
no branch points, and no essential singularities. 

The purpose of this paper is to find the function 
sr(po,p) whose Fourier transform has the asymptotic 
behavior in Eq. p. 1(1 . The physical approximations 
which lead to Eq. I|1.1H are valid in the region po ~ p 
000 El E3 and it is only this region that is inves- 
tigated here. Sec. II shows that the Bloch-Nordsieck 
propagator Sjt(t,p) has the very unusual property of be- 
ing an entire function of complex t. In Sec. Ill only the 
leading term in the asymptotic expansion is employed, 



aT 



■ exp 



■ Po-P 

i 

2aT 



a exp 



■ Po-P 



aT 



(1.4) 



This approximation for sr(pq,p) works surprisingly well. 
It is checked by showing that the Fourier transform over 
energies po ~ p agrees with Eq. (|l.lfl at large t. The mo- 
tivation for the ansatz, Eq. (|1.4|) . is given in Appendix 
B. There a complete asymptotic expansion of sr(pq,p) is 
obtained. The first term of this expansion is Eq. (|1.4J) . 
Sec. IV shows how the absence of singularities in the elec- 
tron propagator imply that in hard scattering processes, 
there will be no infrared or collinear divergences and the 
quantitative contribution of these regions is numerically 
small. Sec. V discusses some further consequences. 



II. COMMENT ON COMPLEX TIME 

The logarithmic dependence on t in Eq. (f L . If) does not 
imply that there is a branch cut for negative t. The full 
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Bloch-Nordsieck calculation of Blaizot and Iancu 
gives 

s R (t,p) = -i6(t)e- ipt exjp[-aTtf(Lo p t)]., 
where f(z) is given by 



ds ■ 



f(z)=C-j + 



This representation shows that f(z) is analytic every- 
where in the complex z plane. (See also Appendix B.) 
Thus the retarded propagator with 9(t) omitted is ana- 
lytic everywhere in the complex t plane. 

Analyticity in t is very unusual. Appendix A shows 
that in general the retarded propagator (without the 0(t)) 
is the difference of two functions: one, S>(t,p), is an- 
alytic in the open strip —/3 < Im (t) < 0; the other 
S < (t,p) is analytic in the open strip < Im(i) < [3. 
Thus generally, the retarded propagator is only defined 
on the real t axis and is not analytic off-axis. Appendix 
A also shows that the analyticity in t of s R (t,p) does 
not improve on, or detract from, the validity of the 
Kubo-Martin-Schwinger condition [12J, which requires 
S > (t-0,p) = -S < (t,p). 



III. THE PROPAGATOR FOR p » p 

This section will provide evidence that a good approx- 
imation to the retarded electron propagator in the region 



P 



naT ixaT 
—— < Repo < p-\ — , 



is the function 
sr(po,p) = -iNexp 



2aT 



aexp 



Po-P 



aT 



(3.1) 
, (3.2) 



where N — \j2-naj aT . It is and entire function of po and 
has the correct behavior anticipated in Eq. I|1.2|l . The 
Fourier transform of Eq. i|3.2|) will be shown to reproduce 
Eq. at large time. The motivation for the ansatz 

(|3.2|) is contained in Appendix B. That appendix contains 
a complete asymptotic expansion of s R (po,p) valid in the 
strip (|3.1|l . Eq. I|3.2|) is first term in the expansion. 
For any propagator, the spectral function (for real po) 

is 

p(Po,p) = i[sii(pa,p) - sr.(po,p)*}- 
For the approximate propagator Eq. I|3.2|l . 



p(po,p) 



2N exp | a cos 
'Po-P , 



Po-P 
2aT 



X cos 



2aT 



Po-P 



aT 



If a is small, this is positive in the range given by Eq. 
(|3.1|) . For |pq — p\ <C airT the approximate behavior is 



p(p ,p)^2Ne a 1 



Po-P 



aT 



A. Context for application 

When the fermions are massless in QCD or QED, chi- 
rality and TCP invariance guarantees that the matrix 
structure of the retarded fermion propagator is 

Sr(po,p) = 7^(7° -7 ■p)sb,{po,p) 

- \(7 +i-p)[sr{-p*o,p)]*- 

In applications, the fermion propagator appears in com- 
bination with a fcrmionic current of the form 



Jx(x')=gTr 



(A^x l ) 1 ^ x {x')0{y 1 ,...y n 



1 

where Z = Tr[exp(— (3H)] is the partition function. Af- 
ter projecting out the Dirac matrices from S a \^J\ a typi- 
cal application the retarded propagator will produce the 
combination 



9(t,p) = -i f dt' s R (t ~ t', P ) J(t', P ). 

J — OO 



(3.3) 



is the response at time t to the source function J acting 
at an earlier time t' . This can be expressed as an energy 
integral 



dpo 
2ir 



e- tpot s R ( Po ,p)J(po,p), (3.4) 



where J(po,p) is the Fourier transform of J(t,p). 

Example: For the conventional retarded propagator 
of the quark type (with a finite 7) the necessary energy 
integration is 



dp 



-ipot 



2tt pq - p + i^f/2 



J{P0iP) 



If J{po,p) is analytic in po in the lower-half plane and 
vanishes sufficiently rapidly as Im pq — > —00, then the 
Po contour can be closed in the lower-half plane and the 
result evaluated by Cauchy's theorem. The simplest way 
to insure that the source has these properties in pq is to 
insist that it vanish for times later than some if. 



J(t',p) = 







t'>t f 



arbitrary t' < tf. 
The Fourier transform of the source is 



J(po,p) 



dt' e m,t 'j(t',p). 



(3.5) 



(3.6) 



This is analytic for Im po < and so the po integration 
to be performed by Cauchy's theorem to give 



t > t 



f 



*(t,p) = exp ( 



ip - 



1\ 



)t J(p- 



i\,P)- (3-7) 



The remainder of this section will show that with a cur- 
rent of this type, the simple function given in Eq. I|3.2|) 
produces the analogous result but with 7/2 replaced by 
aT[ln(uj p t) + C}. 



3 



B. Saddle point integration 

The asymptotic behavior in Eq. 1|1.1|) was obtained in 
0, H, 01 E3 by retaining infrared dominant contributions 
in the vicinity po « p. The ansatz in Eq. i|3.2[l is to 
be tested in the strip l|3.1|) . To isolate the behavior of 
the propagator at po rs p, it is necessary to introduce 
some smearing function f(po) that suppresses the effects 
of large po — p. With the smearing, Eq. I|3.4fl becomes 



dp 
2T 



- ipot s R (po,p)J(po,p)f(po) 



Substituting the ansatz Eq. (|3.2I) gives 
dp 



V(t,p) = -iN 
where 

<P(Po) 



2tt 



e^J(p ,p)f( Po ), (3.8) 



ip t + i 



: Po-P 
2aT 



a exp 



■ Po-P 



aT 



(3.9) 



As t — > oo the integrand oscillates so rapidly that most 
values of po contribute very little. The dominant contri- 
bution comes from the region in which <j)(po) is stationary. 
The first derivative of d> is 



1 



<t>'(p ) = -it + i— + 



; exp 



Po-P 



2aT aT 
The stationary point p satisfying 4>'(p ) 

'aTt - 1/2 



aT . 
= is 



p a = p — iaT In 



(3.10) 



An increasing value of t moves the saddle point further 
down into the lower half of the complex plane. At the 
saddle point, </> itself has the value 





'aTt -1/2" 






a 





0(p o )= -ipt + (aTt - 1/2) ^ - In 



The Taylor series expansion of 4>(pa) about the saddle 
point is 



<t>(pa 



<t>(Po) + 2^"(^o)(Po - Po) 



+ 



The necessary second derivative is 

'aTt -1/2 



4> (pq) = 



The fact that the second derivative is real and negative 
means that the value of Re [</>(po)j is a maximum at the 
saddle point and that the value decreases if po — p is 
positive real or negative real. The dominant contribution 
to the integral is 



tf(t,p) = -iN exp[(f>(p )} 



dpo 
2tt 



J(pa,p)f(Po) 



x exp<j - -(aTt - 1/2) 



Pa 



Po l2 



aT 



The integration contour can be shifted down below the 
real po axis. At large time aTt 3> 1, if the p contour is a 
straight line parallel to the real axis keeping the difference 
Po — Po rea l the integrand falls off very rapidly. The 
integration yields 



*(*,P) 



N 



exp[<Kp )]- 



aT 



^faTt- 1/2 
Using the value of 4>(p ) and N gives 



J(P^P)f(Po)- 



V(t,p) 



-i exp ■ 



1 




raTt - 1/2- 






aTt 


In 


-1 


} 


2 




a 







x J(Po>P)f(Po) 



The relation ln(aT/a) — \n(oj p ) + C + 1 is useful in sim- 
plifying the result at large time, aTt ^> 1/2, and gives 
the final form 



^(t,p) = — i exp | — ipt — aTt 
xJ(Po'P)f(Po)- 



ln[uj p t] + C] } 



(3.11) 



The time dependence coincides with Eq. (|l.l(l . The sta- 
tionary point is p = p — iaT[\n(ujpt) + C + 1] at large 
times. This confirms the ansatz. 



IV. SOFT AND COLLINEAR PHOTONS 

Knowing the electron propagator makes it possible to 
compute radiatiave processes. For a massive electron in a 
QED plasma there are various infrared divergences asso- 
ciated with real and virtual photons that eventually can- 
cel. The divergences arise because if an on-shell electron 
with P 2 — m 2 emits an on-shell photon with K 2 = 0, 
the electron subsequently propagates with amplitude 

1 -1 



(P-K) 2 



2\\s\(E -pcost 



(4.1) 



For very soft photons, |k| — * 0, this can produce an in- 
frared divergence. For example, when an electron of mo- 
mentum P undergoes a hard scattering to momentum 
P' , infrared divergences arise from the real emission and 
real absorption of soft photons by the incoming electron 
and by the outgoing electron. Infrared divergences also 
arise from virtual photons in three ways: from electron 
self-energy corrections to the ingoing electron and to the 
outgoing electron; and from virtual photons that link the 
incoming and outgoing electrons (i.e. corrections to the 
hard scattering vertex). For massive electrons, the in- 
frared divergences from the real photons exactly cancel 
those from the virtual photons |13l Il4j . 

If the electron is massless then E = p in Eq. 1)4. 1|) 
and there would also be a collinear singularity arising 
from 1/(1 — cos 9). As Blaizot and Iancu showed 0, 
the thermal propagator for massless electrons does not 
actually have a pole at po — p; it is instead an entire 
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function of pq. Since there is no pole, there can be no 
infrared divergences and no collinear divergences. 

However, even though there are no divergences there 
is still an important quantitative question of how large 
the radiative corrections are. To compute these correc- 
tions, the first step is to compute the vertex that couples 
a photon to the electron propagator. At this stage it is 
convenient to return to the strict Bloch-Nordsieck struc- 
ture of @, H, in which the propagating electron has 
a velocity vector v. Thus in Eq. 13.2fl the combination 
Po — p should be replaced by p — v ■ p = P ■ v. The Ward- 
Takahashi identity for a photon of four-momentum K to 
be radiated by an electron with initial momentum P and 
final momentum P — K is 

- K^(P -K,P) = a£{P -K)- s£{P). (4.2) 

The general solution for the electromagnetic vertex is 



t^(p-k,p) 



K-v 



(4.3) 



(In the usual case of a free electron s R ^{P) = P ■ v so 
that (P — K,P) — u M as appropriate for the Bloch- 
Nordsieck approximation.) 

Now we employ this in the example of an electron un- 
dergoing a hard scattering from P to P'. The initial and 
final momenta are on shell; P-v — and P' -v' — 0. Dur- 
ing the hard scattering, either the incoming electron or 
the outgoing electron can emit a real photon with four- 
momentum K. If Mb aro (P',P) is the amplitude for the 
hard scattering without radiation, the amplitude for hard 
scattering with radiation is approximately 

M hare (P',P)e fi J^(K). 

The effective current to which the photon couples is a 
combination of the electromagnetic vertex above and an 
off-shell electron propagator: 

J*{K) = s R (P-K)gr»(P-K,P) 

+ gT^(P',P'-K)s R (P'-K). 

Substituting for the vertex functions gives 



K-v 



sr(P-K) 
sr(P) 



gv 



K-v 



sr(P'-K) 
■sr(P') 



(4-4) 

When the momentum transfer Q 2 — (P — P 1 ) 2 is large, 
the cross section for hard scattering with the emission of 
one real photon [ko = k) is 



da 




dQ 2 


J 

bare « 



d 3 k 



(27r) 3 2fc 



(l+n^lv/^)! 2 , (4.5) 

pol 



where n = l/[exp(/3fc) — 1] is the Bose-Einstein function. 

For a free electron propagator s^ 1 (P) = and 
s r 1 (P') = so that the current is 



J»(K) 



-gv' 1 



gv 



This would produce collinear singularities at k ■ v = 1 
and fc • v' = 1 as well as logarithmic and linear infrared 
divergences from J d 3 k/k A . 

However, since the electron propagator is actually an 
entire function, it has no poles. It is a function of the 
variable P-v. For small values of K ■ v the current l|4.4() 
becomes 



K ■ v <€. olT : J»(K) 



v^ dsj 



-fJ 



+9 



sr dp 
v" 1 ds R 



s R dp' Q 



Pa=P'-v' 



(4.6) 



In the emission cross section (|4.5|l there is no collinear 
divergence. The integrand becomes J d 3 k/k 2 so there 
is no infrared divergence. The contribution of soft k to 
multiplicative factor in Eq. (|4.5J) is approximately 



d 3 k T ( g dsn 



(2n) 3 2k k \sr dpo 



U (1 + 2a) 2 



aT 



-in 



k(l-k-v) k(l-k-v') 



Since /c max *C aT this is a small correction. A similar 
estimate applies to the effects produced by the absorption 
of a photon and by the exchange of virtual photons. 



V. DISCUSSION 

The approximate propagator in Eq. (|3.2|l has been 
shown to produce the correct large time behavior in Eq. 

This approximate propagator is the first term of 
the complete expansion of sr(pq,p) that is performed in 
Appendix B. The approximate propagator Eq. (|3.2(l has 
some deficiencies. First, the asymptotic behavior ln(cj p t) 
is only reached at a very large time t satisfying aTt 3> 1, 
whereas the asymptotic form in Eq. i|l.lfl se ts m a t a 
smaller time tu p t 3> 1. Second, the corrections to the 
asymptotic time dependence are of order 1 /aTt whereas 
the corrections to Eq. are of order exp(— ui p t)/u) p t. 

The absence of any singularity in s r (pq,p) eliminates 
infrared and collinear divergences in radiative processes. 
However, there are unsettling consequences in higher or- 
ders. For example, the contribution to the vacuum po- 
larization tensor, (Q), of a virtual electron-positron 
pair would normally have a branch cut that signals the 
production of a real e + e _ pair. If sn(po,p) is entire, 
it appears there can be no e + e~ production. This sit- 
uation arises in zero-temperature models of color con- 
finement. There the quark and gluon propagators are 
thought to have no simple poles and perhaps no singu- 
larities at all. A clear example of this occurs in the non- 
relativistic quark model. There the confining potentials 
used in the Schrocdinger equation produce wavefunctions 
that fall with distance like exp[— /(r)], where at large dis- 
tance f(r) grows more rapidly than r. Form factors and 
other matrix elements computed with these wavefunc- 
tions are entire functions of momentum as a consequence 
of confinement 15]. 
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APPENDIX A: COMPLEX TIME AND THE KMS 
CONDITION 

The Kubo, Martin, Schwinger condition ^| requires 
periodicity in complex time of a certain two point func- 
tion S > (x) defined below. This section shows that the 
KMS condition is essentially unrelated to the behavior 
of the retarded propagator in complex time. 



Definitions 



The basic two-point function is 



(S>(x)) 



a/3 



-m 



where g = exp(— /3i7)/Tr[exp(— f3H)\ is the density oper- 
ator. By inserting a complete set of states between the 
operators, it is easy to show that this is an analytic func- 
tion of time in the open strip —j3 < Im (x°) < 0. The 
related function 



is analytic in the open strip < Im (x ) < j3. The two 
functions are related by the Kubo, Martin, Schwinger 
condition |12| 

lm(x )>0: Syix -tj9,x) = -5<(a: ,x). (Al) 

The retarded and advanced propagators are related to 
the thermal average of the anticommutators: 



S R (x) = Oix^S^x) - S^x) 
S A (x) = 9{-x°)\ -S > {x)+S < (x) 



The only common region of definition for S > (x) and 
S<(x) is for x° real. The retarded and advanced propa- 
gators, with (9(±x°) omitted, are usually restricted to x° 
real. 



2. Momentum space 

For massless fermions, the combination of chirality and 
TCP invariance makes the propagators a linear combina- 
tion of 7° and 7 • p: 



MP) = ^-i-p)s* r {pI,p) 

- 7^7° + 7-P)sr(-Po,p). 

The propagators are thus determined by the one function 
sr(po,p). Knowing this, one can also construct the two- 
point functions using 



MP) 
MP) 



1 



1 + exp(-#p ) 
-1 

1 + exp(Ppo) 



S r {P)~Sa(P)\ (A2) 
S r (P)-Sa(P)]. (A3) 



Since the retarded propagator is analytic in the upper- 
half of the complex po plane and the advanced is analytic 
in the lower- half, generally S > (P) and S < (P) are not 
analytic anywhere. Their arguments, po, should always 
be considered as real. 

The Fourier transforms of Eq. |jX2l and auto- 
matically satisfy the KMS condition ilAlfl . The condi- 
tion imposes no restriction on the retarded propagator 
in momentum space, Sr(P). 



APPENDIX B: COMPUTATION OF s R (p ,p) 

This section computes the Fourier transform of Eq. 
and obtains a complete asymptotic expansion of 
sr(P0iP) valid in the strip 



naT iraT 
—— < Re(p )-p< ——. 



(Bl) 



The final result is given in Eqs. (|B9|I and (|B12(1 . 

Before beginning this analysis it is useful to examine 
the full Blaizot and Iancu 0, H, 111 result for the retarded 
propagator: 



sr(po,p) = ~i / dte w ^ 



W(p ,t) = i{ Po -p)t-aTtf{u p t) 



/(z)= C -7 + 



ds 

s 



(B2) 
(B3) 
(B4) 



The integral that defines f(z) can be computed as a 
power series: 

««'- c -' +1 -£l^i5i- 

This series converges everywhere in the complex z plane 
and thus f(z) is an entire function. In particular, f(z) 
has no branch cuts. However the power series expansion 
is useless for large z. One can rewrite Eq. (|B4|) using 
repeated integration by parts to obtain 

1 e~ z (-V\ k W 

m = H z)+c+- +—y / lai± 



z z * — ' z 

fc=i 



:(7° ' + 7 ' P) s*r(-P*o,p) 



+ (-l)" +i (n+l)! 



ds ■ 



-,n+2 ' 



This is also an exact result. Despite the explicit ln(z), 
there is no branch point: If z is moved in the complex 
plane along a counter-clockwise path that encircles the 
origin, then ln(z) changes by 27r? and the integral over 
s changes by —2iri. Thus j{e 2m z) = f(z): there is no 
branch cut. For Rez > and large, one can use the 
approximate form 

f(z) » ln(z) + C. 

As observed by Blaizot and Iancu § , the integrand in 
Eq. l|B2jl falls so rapidly at large t (large z) that the inte- 
gral converges any complex value of po. Furthermore all 
the derivatives d n SR(po,p)/dpQ are finite at any complex 
value of pq. Thus sr(pq,p) is an entire function of pq. 



1. Exact computation 

With the simplified integrand, the Fourier transform 
to be computed is 

poo 

8 R (po,p) = -i dte w{ ^\ (B5) 
Jo 

where 

W( Po , t) = ^(p - p)t - aTt[\n{u p t) + C] . (B6) 

The computation of the propagator at large t in 0, H, 
0, is most reliable at po ~ p. This will be the range 
in which the Fourier transform is performed. It will be 
necessary to require that po satisfy Eq. i|Bl(l . 

The integration contour in Eq. 1|B5|) can be distorted 
into the complex t plane and the convergence will not be 
compromised provided that Re t — > +oo at the terminus 
of the contour. Define the complex time 



-C-i 



t = 



exp 



Po-P 



aT 



(B7) 



In the complex t plane, the time derivative of W vanishes 
at t. The integration contour in Eq. I)B5|1 can be taken 
as a straight line from the origin through t and since 
Ret > the integration converges. The straight line is 
parametrized by a real variable u such that t = tu. The 
exponent becomes 



W(p ,t) = W[-uhx u + 



where 



W = a exp 



.Po-P 



aT 



(B8) 



and a is given by Eq. It will be convenient to 

express W as 

W(p ,t) = W -W[ulnu-u+l] 

The quantity in square brackets vanishes at u = 1 and 
is otherwise positive throughout the domain < u < oo. 



From u = to u = 1 it decreases monotonically; from 
u = 1 to u = oo it increases monotonically. 

Using dt ~tdu and t — W /aT, the retarded propaga- 
tor can be expressed as 



sr(po,p) 



-iW w 



aT 



{l + j}, 



(B9) 



where the integral has been split into two parts: 



du e 



-W[u\nu— u+1] 







J = J due~ w[ulnu - u+1] . 

In the first integral change variables from u to x where 

0<it<l: -x 2 = u < \n(u < ) - U< + 1. (BIO) 

As u increases from to 1, x decreases from y2 to 0. 
Expand u in a series 

oo 

u< =l + J2(-l) n c n x n . 
n=l 

The first few coefficients are Ci=l, c 2 =l/6, c 3 ~— 1/72, 
C4= 1/270. The coefficients satisfy the recursion relation 



n 

l -2 



n+l 



— ^ k Ck Cn 



+2-k- 



k=l 



Then / is easily computed 



I = - 



V2 



du 



dx ^< e -w x </2 



dx 



n=l 



V~2 



dxx^e-™* 2 / 2 . 



In the second integral change variables to y where 

1 < u < oo : iy 2 = u> ln(u>) - m> + 1. (Bll) 

As u increases from 1 to oo, y increases from to oo. 
Expand u in a series 

oo 

71=1 

with the same c„. Then J is 
du 



J 



dy rr> e -«V/2 
o dy 



E 



nc rl 



n-l-Wy 2 /2 



dyy n ~ L e 
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In the sum I + J there is a partial cancellation. The 
integrations can be performed exactly with the result 



J = 



£ (2fc+1)C2fe+i (- 2 dl) {w 



oo f d \ k 

+ E(2 fc + 2) C2fe+2 (-2-=j _ r . 

k— 



erfc(VW^ 
(B12) 




The asymptotic behavior of the complementary error 
function is 



erfc 



Substitution of Eq. i|B12(l into Eq. I|B9|) gives a complete 
asymptotic expansion for the retarded propagator when 
Po satisfies (IB1|I . Although this expansion is periodic 
under po ~ > Pa + 47raT, such a shift is well outside the 
domain of validity given in Eq. I|B1|) . 



2. Large W limit 



The asymptotic series (|B12|) obviously behaves very 
badly at small values of W . The most useful applica- 

tion of the results is clearly when Re W is large and 
positive. In this region the exponentials are negligible 



s R {po,p) = -i 



V2nW 
aT 



E 

k=0 



(2fc + l)!! 

C2k+1 Z=k l B13 J 



-« — — e 



24 

w 



w 

23 
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The simplest approximation is to keep only the first 
term in Eq. (|B13|) . In terms of po and p the propagator 



sr(po,p) 



aT 



■ exp 



Po-P 
2aT 



a exp 



Po-P 



aT 



(B14) 

This is analytic for any finite po ; it vanishes in the limit 
Im po - * +oo; and it diverges when Repo = P and 



Im po — ► — °o in agreement with Eq. (|1.2|) . This is the 
form used throughout the paper. 

Before concluding, it is necessary to discuss how W 
can be large. From Eq. 1B8(I the absolute value is 



W = ae 



-Impo/ctT 



3£ 
4tt 



exp 



C- 1 



Impo 



(B15) 



For small g there are two ways in which the magnitude 
\W\ can be large. One possibility is to restrict consider- 
ation to the region in which Impo/aT is large and nega- 
tive. This is perfectly valid, but it excludes the physically 
important region of real po- The other possibility is for 
the constant C to be large and negative. 

Both groups [U showed that value of C in the 
Coulomb gauge is C = In V3 + 7^ — 1 = 0.12652 ... In co- 
variant gauges the value of C can be different. When the 
zero temperature, Feynman propagator for the photon is 



1 



<r + (i-a) 



it is necessary to impose an infrared cutoff /i. Blaizot 
and Iancu |3] discuss two limits: (j,t 1 and /it ^> 1. If 
fit 3> 1 then C is independent of the gauge parameter 
A. However in the calculation of Sec. Bl, the dominant 
contribution comes from finite values of t. (The integral 
/ sums t in the range < t < t. If the integral lim- 
its on J are changed to < y < y/2, corresponding to 
t < t < et, the error made is 2.35 x e~ w /W, which is 
negligible.) Since the important values of t are finite in 
the previous integration, then as the infrared cutoff [i is 
reduced the appropriate limit is fit -C 1. In this limit 
Blaizot and Iancu [3| showed that C does depend on the 
gauge parameter A: 



0.12652. 



(B16) 



Consequently a large, negative A will make C large and 
negative, which will make \W\ large thus justifying the 
approximation. In [b|| the range —857 < A < 433 was 
explored. 
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